Qubit-assisted thermometry of a quantum harmonic oscillator 



Matteo Brunelli, 1 Stefano Olivares, 1,2 Mauro Paternostro, 3 and Matteo G. A. Paris 12 

1 Dipartimento di Fisica, Universitd degli Studi di Milano, 1-20133 Milano, Italy 
2 CNISM, UdR Milano, 1-20133 Milano, Italy 
3 Centre for Theoretical Atomic; Molecular and Optical Physics, 
School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 INN, United Kingdom 

(Dated: March 4, 2013) 

We use the theory of quantum estimation in two different qubit-boson coupling models to demonstrate that 
the temperature of a quantum harmonic oscillator can be estimated with high precision by quantum-limited 
measurements on the qubit. The two models that we address embody situations of current physical interest due 
to their connection with ongoing experimental efforts on the control of mesoscopic dynamics. We show that 
population measurements performed over the qubit probe are near optimal for a broad range of temperatures of 
the harmonic oscillator. 
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I. INTRODUCTION 

The improved control over systems of intrinsic complex- 
ity makes the implementation of techniques for the inference 
of specific properties of their states a necessary step towards 
the achievement of full quantum control. Yet, it is often the 
case that the device into which we would like to enforce quan- 
tum mechanical features is not fully or easily addressable. 
System-interrogation can only be performed, in such cases, 
in an indirect way through the use of probes of an appropri- 
ate nature HI. It is thus very important to devise experimen- 
tally implementable strategies for the inference of properties 
of inaccessible quantum systems, identify the optimal state- 
preparation of the probe as well as the observable that allows 
for the maximum extraction of information about the parame- 
ter that we are interested in. 

This agenda is dressed of even more relevance due to the 
recent experimental efforts produced towards the quantum- 
limited management of mesoscopic systems, such as super- 
conducting devices [2 J, light-interfaced cold-atom systems [3 1 
and mechanical systems operating at the quantum level Q. 
All such systems have in common the use of "quantum in- 
terfaces" with devices of a different nature, which are then 
exploited for state-preparation, manipulation of information 
and possibly read-out. The quantum-interface paradigm is 
indeed very fruitful for the extraction of information out of 
a system that is only partially accessible: through the cou- 
pling with a controllable subsystem, one can indeed arrange 
for mechanisms able to provide useful knowledge on key fea- 
tures of a dynamics or a state. Examples of such a possibil- 
ity, which have been materialised in successful experimental 
demonstrations, include the micro-maser technology for the 
revelation of the properties of the field within a high-quality 
microwave cavity [5 |, the coupling of a Bose-Einstein con- 
densate to a (classical) mechanical oscillator for the investiga- 
tion on Casimir-Polder effects [6|, and intra-cavity quantum 
optomechanics, where the radiation-pressure force is used to 
read the noise properties of a mechanical mode [7]. In the 
latter context, in particular, a key parameter is embodied by 
the temperature at which the mechanical mode operates. In- 
deed, unwanted thermal effects typically spoil the quantum 



features enforces in the mechanical system by means of a co- 
herent evolution. Having a precise quantitative estimate of the 
entity of such effects |8| would be crucial not only for pre- 
diction purposes but also to design in the best possible way a 
quantum-enforcing protocol that accounts, ab initio, such un- 
desired effects. Needless to say, these considerations can be 
extended straightforwardly to any of the scenarios addressed 
above. 

Recently, strategies for the determination of the tempera- 
ture of a harmonic oscillator have been put forward, based 
on the coupling to a quantum probe embodied by a two-level 
system (a qubit) (9) \W\. The coupling model to be sued 
for the thermometry of the oscillator's state was the Jaynes- 
Cummings one, within and beyond the so-called rotating wave 
approximation ifTTl . While Ref. ifTOll proposed the use of the 
ac Stark effect as a way to infer the temperature of the oscilla- 
tor, Brunelli etal. [ 9 ] have applied the proper tools of quantum 
estimation theory (QET) ifTZHHl to design optimal protocols 
for the estimate of the system's temperature. 

In this paper, we significantly extend the approach in |9| 
to other physically motivated qubit-oscillator models, prov- 
ing that optimal and effective thermometry can indeed be per- 
formed by means of simple measurements onto the qubit's 
state. We tackle both the coupling between a superconducting 
qubit and a nano electromechanical oscillator and the far-off 
resonant interaction between a two-level atom and the field of 
a cavity, thus providing an analytical QET-based study of an 
ample spectrum of experimentally motivated situations. 

The remainder of this work is organized as follows: in 
Sec. [II] we describe the general system that we address and 
introduce the QET tools for our analysis. Sec. [Til] studies 
the first model of our investigation, which addresses the ca- 
pacitive coupling of a superconducting qubit and a nano- 
electromechanical oscillator. In Sec. [lV]we assess our QET- 
based approach in the case of a qubit that is off-resonantly 
coupled to a harmonic oscillator, such as for a two-level atom 
in a far-off resonant cavity. In both instances, population mea- 
surements over the probing qubit allow for the optimized es- 
timate of the oscillator's temperature. Finally, in Sec. [V] we 
summarize our findings and open up new perspectives. 
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II. APPROACH TO THE PROBLEM 

Let us consider a general quantum harmonic oscillator with 
frequency and at thermal equilibrium with its environment. 
The state of the oscillator is described by the Gibbs density 
operator (we use natural units, i.e. H = 1, throughout the 
manuscript) 
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where n = (e^ n — 1) 1 is the average number of thermal ex- 
citations, |n) is a Fock state with n quanta, Z = Tr[e~^ na a ] 
is the partition function and a (a}) is the bosonic annihilation 
(creation) operator of the harmonic oscillator. Our aim is to 
estimate the inverse temperature /3 = 1/fceT of the oscilla- 
tor by coupling it with a qubit encoded in the logical states 
{\0) q , \ l) q } of a two-level system that is initially prepared in 
the general pure state 
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Here (0, ip) are the angles defining the orientation of the 
qubit's Bloch vector in the corresponding Bloch sphere, while 
ks is the Boltzmann constant. We assume no initial correla- 
tion between the probe and the oscillator and also assume that 
the interaction Hamiltonian has the general form 



Hj = g A Q <E) A 



(3) 



where A Q (A q ) is an operator in the Hilbert space of the oscil- 
lator (qubit) and g a coupling constant. In what follows, we 
shall call {|x) } a basis of states of the harmonic oscillators 
that are eigenstates of A Q , i.e. A Q \x) a — x \x) Q . 

Any measurement aimed at estimating the temperature of 
the oscillator is performed on the state g q of the probing qubit 
after its joint evolution with system o. That is 
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where p (x) = (x\p \x) are the diagonal matrix elements 
of the initial thermal state in the basis |a;) of the oscillator 
operator A Q . In what follows, we make use of the apparatus 
of QET to design the optimal probing state and measurements 
needed to estimate the inverse temperature (3. According to 
the Cramer-Rao inequality, the variance 5(j) of any unbiased 
estimator of an arbitrary quantity 7 satisfies the inequality 
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with M the number of measurements used in order to perform 
the estimate and F(j3) the Fisher information of /3, which is 
defined as 



Pi 



(5) 



where pj is the probabilities to get outcome j from a measure- 
ments performed over the qubit state and described, in gen- 
eral, by the positive operator valued measurement (POVM) 
{IL, : IT, > 0, J^jiij = 1}. Such probabilities are calcu- 
lated assuming the oscillator at the inverse temperature f3, i.e. 
Pj = TrMfilLj}. 

The quantum mechanical counterpart of the Fisher Infor- 
mation is defined as 



H(f3)=Tr{Q q L 2 (!3)} 



(6) 



with L(j3) the symmetric logarithmic derivative operator, sat- 
isfying the equation 



dpg q = [L(p) Qq + g q L((3)}/2. 



(7) 



The quantum Fisher Information (QFI) is an upper bound for 
F(j3) as it embodies the optimization of the Fisher Informa- 
tion over any possible measurement performed over the prob- 
ing qubit states. The QFI is thus independent of the specific 
measurement strategy and is an intrinsic feature of the family 
of probing states. Eq. Q can then be rewritten as 

which extends the Cramer-Rao bound to the quantum domain 
and embodies the ultimate limit to the precision of the esti- 
mate of p. A measurement is optimal when the correspond- 
ing Fisher information F(f3) equals the quantum Fisher Infor- 
mation H(j3). Although various instances of optimal mea- 
surement may be found, depending on the model at hand, 
the observable embodied by the spectral measure of L(j3) is 
certainly optimal. Upon diagonalization of the probe state 
Q q = 9+ (V>+I g + Q- IV 1 -) (ip-\ r the QFI can be com- 
puted explicitly as 
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with 7 = (1 — 2g+) 2 [?)]. In what follows, we consider two 
exactly solvable models corresponding to interesting phys- 
ical situations and compute the QFI to assess the ultimate 
precision in the estimation of temperature achievable by any 
measurements performed on the states of the probing sys- 
tem. We will compare such optimal performance to what 
is obtained through the Fisher Information associated with 
population measurements of the probe, i.e. for {IT,} = 
{|0) (0| , |1) We show that in some cases, population 

measurements are indeed optimal for the estimation of tem- 
perature. 



III. JAYNES-CUMMINGS COUPLING BEYOND THE 
ROTATING WAVE APPROXIMATION 

The first model that we address corresponds to the choice 
A a = X a = (a + a')/y2, i.e. the in-phase quadrature op- 
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erator of the harmonic oscillator, and A q — a x , which is the 
x-Pauli spin operator. Correspondingly, the interaction reads 



H, 



This model is encountered in a few different contexts. On 
one hand, it describes the effective interaction Hamiltonian 
for the electric-dipole coupling between a two-level atom and 
the field of a cavity, thus embodying the celebrated Jaynes- 
Cummings Hamiltonian [11 J beyond the so-called rotating 
wave approximation. Moreover and rather less intuitively, 
the same model is achieved by considering a nanomechani- 
cal oscillator (a nano beam) coupled capacitively to a Cooper- 
pair box (CPB) operating at the so-called charge degeneracy 
point 1161 . where the dynamics of the CPB can be righteously 
be approximated to that of a two-level system encoded in the 
space spanned by states | ±) = (|0)±|l))/%/2. Here {|0) , |1)} 
are states with exactly and 1 excess Cooper pairs in the large 
superconducting island shown in Fig. [T] The natural Hamilto- 
nian of the system reads 



Hi = 
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such that Q g ~ 0, the interaction Hamiltonian of the system 
can be cast into the form 



Hi = X(a 
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with A = g/V2 an effective coupling rate whose form is 
inessential for our tasks. 

The estimate of the temperature in this particular context 
is especially relevant. Indeed, the nano beam is in contact 
with a thermal phononic background due to the substrate onto 
which it is nano-fabricated ifTTl . The coupling with the su- 
perconducting qubit addressed above holds the potential to 
prepare non-classical states of the nano-beam. Indeed, the 
time-evolution operator corresponding to Eq. (Ill reads, in 
the qubit basis, as 

U(t) = cos(gtX)i q — ism(gtX)a x . 



Let us assume that the nano-beam is initialized in a coherent 
state |a) (a € R), while the qubit is prepared in \+) q . The 
evolution will generate the qubit-oscillator state 
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\ri(a)) qo = ^(\a + igt) o \0) i 
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with Q and <p the canonical charge and phase operator of the 
CPB, Ct the total capacitance of the island, Q g — C g V g + 
C X V X the total gate charge, Ej the Josephson energy and 
the frequency of the nanomechanical oscillator, as indicated 
above 061 . By defining a x = |+) (— | + |— ) (+|, expand- 
ing Hi in series of the ratio x/d between the actual position 
of the mechanical oscillator and its equilibrium distance from 
the CPB (the amplitude of the oscillations is assumed to be 
small enough that only terms proportional to x/d are retained 
in such expansion) and adjusting the gate and driving voltages 



for gt ~ 7r the two co- 




FIG. 1: (Color online) An electrically driven nanomechanical oscilla- 
tor (bias voltage V x ) is coupled to a CPB through the capacitance C x - 
The state of the CPB is controlled by the gate voltage V g (coupled to 
the box through the capacitance C g ) and the Josephson energy Ej. 
We work at the charge degeneracy point. 



As | (a - igt\a + igt)\ 2 = e~ 4s * 
herent states \a ± igt) are quasi-orthogonal and Eq. ( 12 1 is al- 
most maximally entangled. By projecting the qubit onto \+) q , 
we achieve the coherent-state superposition Af(\a + igt) + 
\a — igt) Q ) (N is a normalization factor) which embodies, in 
the limit of quasi-orthogonal coherent states mentioned above, 
a highly non-classical state. However, a thermal-state prepa- 
ration of the harmonic oscillator will smear out such non- 
classicality, pushing the state towards the statistical mixture 



Po,th= d 2 aG(a,V)\r](a)} 00 (r](a)\ 



(13) 



where G(a, V) is a Gaussian distribution of width V — 
2n + 1 ifTSl . Determining the exact initial temperature of 
the nano beam is thus key for the success of such conditional 
strategies for the enforcement of non-classical features. Our 
approach to the estimate of f3 will follow the general strat- 
egy described above, which we now describe for the specific 
model in Eq. 
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The elements of the state of the probing qubit after the in- 
teraction with the oscillator and the trace over its degrees of 
freedom can be calculated explicitly as 



0,03)= 



1 + cos 6e t sin ( 
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1 — cos 6e~^ 



(14) 

with C = coth {(\ t 2 , where j3 — and t = gt are 

respectively the dimensionless inverse temperature and inter- 
action time. The Fisher Information associated with a mea- 
surement of the populations of g q (j3), i.e. a measurement of 
the z-Pauli operator a z , reads 



cos 2 6>csch 4 (13/2) t 4 
e 2C - - cos 2 6 4 
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which is a function of (3, 9, r. Compared to the case where the 
rotating wave approximation is invoked [9], i.e. for a qubit- 
oscillator interaction of the form g(a)<j_ + h.c.) with a± the 
ladder operators of the qubit, the Fisher Information displays 
a symmetric behavior with respect to 9 and is no longer a pe- 
riodic function of the time r. The maximum is achieved by 
choosing 8 = {0, ir}, i.e. by preparing the qubit either in |0) 
or |1), while for 6 = | the Fisher Information identically 
vanishes. Finally, upon choosing one of the optimal qubit 
preparations, we found that the maximum value of the Fisher 
Information only depends on temperature. In fact, the value 
of the dimensionless time at which the Fisher Information is 
maximized is the following function of /3 
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with W(y) the Lambert function of argument y 112011 . The be- 
havior of T opt {(i) is shown in the left panel of Fig. [5] while 
in the right one we give examples of the form of the F(f3) for 
large temperatures. Remarkably, in the relevant regime ^ > 1 
(low temperatures), the optimal interaction time becomes al- 
most independent on /3. This means that no fine tuning of the 
interaction time is needed and only a rough a priori informa- 
tion is needed to implement the optimal measurement. 

In order to evaluate the quantum Fisher Information, we 
have diagonalized the state of the probe, as described in 
Sec. [II] The explicit calculation, which produces expressions 
too involved to be reported here, shows that H(f3) is maxi- 
mized for two independent sets of choices of the qubit-state 
parameters. One can either prepare the qubit in one of the ba- 
sis states |0) or |1), independently of the angle ip, or choose 
ip = {tt/2, 37r/2}, regardless of 9 (cf. Fig. [3]). The values of 
the QFI are the same in both cases, and the analytic expres- 
sion of H(f3) reduces to the one taken by F((3) for the choice 
9 = 0. This demonstrates that population measurements are 
optimal for the whole range of temperatures. 

This conclusion is further strengthened by the analysis of 
the spectral measure of the symmetric logarithmic derivative 
L(f3). For both the optimal probe-state preparations, L(J3) is 




FIG. 3: (Color online) Temporal evolution of the quantum Fisher In- 
formation for P = 10, ip = and 8 = (blue), 6 = j (magenta), 
8 = |+0.1 (green). For 8 = 0, the quantum Fisher Information 
is maximum (independently on ip) and equals the Fisher Information 
associated to population measurements, while for 8 = ~ H (f3) iden- 
tically vanishes. Inset: the quantum Fisher Information evaluated at 
the optimal time, plotted against 8 and for if = (blue), ip = J 
(red), p = \ (orange). The last choice leads to the same maximum 
value of the quantum Fisher Information and does not depend on 8. 



diagonal and reads 



csch 



(cotK -1)1+ cschC<r 2 



(16) 

with a z the z-Pauli operator. The explicit presence of such 



operator in Eq. ( 16 1 demonstrates the optimality of population 



measurements for the estimation of temperature in this model. 

We end the section by noticing that if the harmonic os- 
cillator is moved away from the equilibrium position, i.e. 
its initial state is described by the displaced thermal state 
D(a)g D> (g), D(a) = expjaa^ — aa} being the displace- 
ment operator, then the probe qubit after the interaction is 
given by 



Q q {a,P) = e 



where g q {(3) is the probe state of Eq. ( 14 1. The QFI is equal 
to the zero displacement case, whereas the Fisher information 
is in general smaller than for zero displacement. 
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FIG. 2: (Color online) Left: The functional form of T opt (/3). Right: 
F(f3 = 3) against 8 and r. Clearly, the Fisher Information is opti- 
mized at 8 = 0, n and quickly decays as r grows. 



IV. FAR OFF-RESONANT SPIN-BOSON INTERACTION 

We now address a second qubit-oscillator interaction 
model, specified by taking A a = a) a and A q = a x . The 
interaction Hamiltonian thus becomes 

Hi = ga'a®a x . (17) 

This model describes for a two-level system interacting far 
off-resonantly with a bosonic mode. Let us consider a two- 
level system (bosonic mode) with transition frequency w (f2), 
interacting through a Jaynes-Cummings model with strength 
A. We call A = 57 — lo the detuning between the two systems. 
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The corresponding time-evolution operator can be written, in 
the basis {|1) , |0)} of the two-level system, as ifTSIl 
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With Q q ,ia = e* q fll, Qq,U = 1 - QqfiO 



and 
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where 



K n = 



sin(f2 n i) 



Q n = yj jA 2 + X 2 a^a is the effective Rabi frequency op- 
erator, and a (a 1 *) are the field operators of the boson. For 



A 2 /4 > X 
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we have 
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and by moving to a reference frame rotating at frequency A, 
we gain the effective picture |[T9l 
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We can now shift the energy of the two-level system so that 
\g) becomes the zero-energy state and 

fj = e~ iA ' 2 * = e -i^«at|l)<l|t 

with H' 2 the interaction Hamiltonian between the boson and 
the two-level system. By reminding that |1) (1| = 2a z + i and 
neglecting an inessential term depending only on the qubits' 
degrees of freedom, we gather the non-trivial interaction term 
2 ^-a 1 'a (g> a z . This is locally equivalent (via a Hadamard gate 
applied to the two-level system) to the model H2 in Eq. ( flT) . 
We now assume this interaction model for the probe-oscillator 
dynamics and the protocol for the estimate of the temperature. 

The matrix elements of the probe state after evolution and 
the trace over the oscillator are 



Qq.oo= cos 2 - + r[2C cos sinc 2 r— sin 9 sin <p sin(2r)], 

sin^cosw n( . o . . . . , ., 

Qq,0i= Hi { [e — cos(2r)Jsm 6 sin i^cos 6 sm(2rj) 
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The expression taken by the Fisher Information is, in this 
case, too lengthy to be reported. Qualitatively, F(j3) depends 
on both 9 and ip and, as in the Jaynes-Cummings model un- 
der the rotating-wave approximation, is a periodic function of 
time r. The probe state preparation that optimises the Fisher 
Information is again 9 = {0, it}. For both such choices, F((3) 
is independent of cp. However, as soon as the qubit initial- 
ization deviates from the optimal cases, F(j3) suddenly drops 
by several orders of magnitude, as shown in the left panel of 
Fig. |4j displaying a weak dependence on ip [taking ip = 5 
gives the maximum of F{(3), see the right panel of Fig. Bf. 
However, the values attained by the Fisher Information in at 
such optimal values of ip are negligible with respect to those 
associated to 9 = {0, tt} [there is a difference of four orders of 
magnitude between the values in the left and the right panels 
of Fig.|4), which makes 9 the only effective qubit parameter. 
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FIG. 4: (Color online) Left: FI for /3 = 10, tp 
the effective time r, for different 8 values: — (blue), 6 — 0.01 
(magenta) and 6 = 0.1 (green). For 9 — the Fisher Information 
does not depend on (p. Right: Fisher Information for f3 — 10 and 
6 = j as a function of time for different <p: ip — | (blue), tp = j 
(magenta) and tp — ^ (green). Although the choice <p — ^ maxi- 
mizes the Fisher Information, it is evident that the relevant parameter 
in setting the qubit is 9. 



If we prepare the qubit in an eigenstate of a z , the Fisher 
Information associated to a population measurement reads 



2e 2/3 sin 2 t [1 + sinh/3 - cos(2r)e~ /3 ] 2 tanh (§) 



/3>i 



F(/3) opt = ^ s £ii e~>" sin 2 t (22) 

(C 3 - 1) [1 + (e/ 3 - cos(2r)) (2e^ - 1) - cos(2r)] [cos(2r) - cosh/3] 2 



On the other hand, by inspecting H(/3) we found, as before, 9 — {0, tt} (regardless of ip), or tp = |7r} (regardless of 
a symmetric behavior with respect to the qubit parameters: at 9) and the values achieved by the QFI are equal in both cases, 
a given dimensionless time r, H(f3) is maximum either for Hence the effective contribution to the dynamics comes from 
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one octant of the Bloch sphere. The analytic expression of the 
optimal quantum Fisher Information is 



reads as follows 



H op t(P) 



sin 2 t [2 cos(2r) - 2 cosh/3 - sinh 2 /3] 



2[cos(2r) - cosh/3] 
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The crucial point here is that, for this model, the optimal 
Fisher Information and its quantum mechanical counterpart 
are no longer the same, i.e. population measurements are 
not the optimal one for the whole set of parameters: at low 
temperatures, provided that we choose an optimal qubit-state 
preparation, we retrieve the optimality of a z measurements as 
Fopt(P) = H opt ((3). When the temperature is raised, on the 
other hand, small discrepancies appear between the temporal 
behavior of these quantities, suggesting that a z is not longer 
the best measurement strategy. This is shown in Fig. [5] where 
we study the relative difference between optimized Fisher and 
quantum Fisher Information. 




FIG. 5: (Color online) Logarithmic plot of the quantity 

{H opt (f3)-F opt (f3)}/H opt (f3) for /3 = 10 (blue), /3 = 5 (magenta) 
and P — 1 (green). 
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(25) 

and shows explicitly that optimality of population measure- 
ment is recovered for in the low temperature regime j3 1. 



V. CONCLUSIONS 

We have addressed the thermometry of a (directly inacces- 
sible) quantum harmonic oscillator through its coupling to a 
quantum probe embodied by a controllable qubit that can be 
subjected to any measurement. By focusing our attention on 
two models of current physical relevance and using the frame- 
work of the (quantum) estimation theory, we have determined 
the preparation of the probe qubit, the measurement and the 
value of the interaction time that optimize the estimate of the 
oscillator's temperature. We found that population measure- 
ments performed over the probing system are nearly optimal 
for an ample range of temperatures. This quite important from 
the operational point of view, given the handiness of imple- 
menting <t 2 measurements in all of the settings that have been 
explicitly addressed here. Our work thus aims at proposing an 
experimentally viable pathway towards the quantum-limited 
inference of the properties of inaccessible quantum systems, 
demonstrating that the paradigm of the coupling with a (fully 
controllable) low-dimensional quantum system is indeed ef- 
fective. We are working towards the extension of this frame- 
work to explicitly open-system dynamics and the characteri- 
zation of the environmental properties affecting the dynamics 
of the harmonic oscillator. 



Some insight comes from the analysis of the symmetric log- 
arithmic derivative, which reads 



L{fi) = L t - L x a z + L z a z 



(24) 



showing the presence of a contribution coming from a term 
proportional to & x , which is responsible for the differences 
between the two estimators. The expression of the coefficients 
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